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Abstract 



Cross sections for elastic and proton-dissociative photoproduction of J/ijj mesons are 
measured with the HI detector in positron-proton collisions at HERA. The data were col- 
lected at ep centre-of-mass energies y^ s» 318 GeV and yfs « 225 GeV, corresponding to 
integrated luminosities of C = 130 pb -1 and C = 10.8pb _1 , respectively. The cross sec- 
J> tions are measured as a function of the photon-proton centre-of-mass energy in the range 

25 < W 1V < 110 GeV. Differential cross sections d<r/di, where t is the squared four- 
momentum transfer at the proton vertex, are measured in the range |t| < 1.2 GeV 2 for the 
elastic process and \t\ < 8 GeV 2 for proton dissociation. The results are compared to other 
measurements. The W 7P and t-dependences are parametrised using phenomenological fits. 
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1 Introduction 

This paper reports a measurement of diffractive J/ip photoproduction in positron-proton inter- 
actions at HERA, ep — \ e J/ip X. For the elastic regime X denotes a proton, whereas for the 
proton-dissociative regime X denotes a proton-dissociative system Y of mass m p < My < 
10 GeV, as depicted in figures [T] 





Figure 1: Diffractive J/ip meson production in electron proton collisions: a) elastic J/ip pro- 
duction in which the proton stays intact and b) proton-dissociative J/ip production in which the 
proton dissociates to a low mass excited state with mass M Y > m p . 



Diffractive vector meson production is characterised by the t-channel exchange of a colour- 
less object between the incoming photon and proton. In the high-energy limit Regge theory 
predicts [1,2] an approximate cross section dependence a oc WjL as a function of the photon- 
proton centre-of-mass energy W 1P . For elastic production of light vector mesons (p, to, (p) 
exponents S ~ 0.22 J3j are observed. In contrast, the cross section for elastic J/ip production, 
IP ~^ J/4>P, ri ses more steeply with W 1P , 5^0.7 PlBl, and is thus incompatible with a univer- 



sal pomeron hypothesis |2J. The W 1P dependence of proton-dissociative J/ip production |6j48| 
is expected to be similar to the elastic case. 

Due to the presence of a hard scale, the mass of the J/ip meson, calculations in perturbative 
Quantum-Chromo-Dynamics (QCD) are possible. The diffractive production of vector mesons 
can then be described in the proton rest frame by a process in which the photon fluctuates into 
a qq pair (or colour-dipole) at a long distance from the proton target. The qq pair interacts with 
the proton via a colour-singlet exchange, which in lowest order QCD is realised as a colourless 
gluon pair [|9- 12 1. The steep rise of the cross section with W 1P is then related to the rise of the 
square of the gluon density towards low values of Bjorken x [|T3|-[T7]]. 



The elastic and proton-dissociative J/ip cross sections as functions of the squared four- 
momentum transfer t at the proton vertex show a fast fall with increasing \t\ [4, 5J [T8| - |25| . For 
the elastic J/ip cross section the ^-dependence can be parametrised by an exponential function 
da/dt ex e~ bel ^ as expected from diffractive scattering. In an optical model the t-dependence of 
the elastic cross section carries information on the transverse size of the interaction region. The 
proton-dissociative cross section falls less steeply than the elastic one and becomes dominant at 
\t\ > 1 GeV 2 . The differential proton-dissociative cross section as a function oft is parametrised 



with a power-law function da/dt oc (1 + (b pd /n)\t\) n , which for low \t\ has an approximate 
exponential behaviour, oc e~ bpd ''L 

Diffractive J/ip production has been studied previously at HERA at low values of \t\ PJHJ 



18-23|, and also at very large values of \t\ [24 25 1, where proton-dissociative J/ip production 



dominates. 

In this analysis cross sections are determined simultaneously for the elastic and proton- 
dissociative regimes. In addition to a measurement at the nominal ep centre-of-mass energy 
of ^fs 7H 318 GeV, data recorded at a lower centre-of-mass energy of y/s ~ 225 GeV are 
analysed. This low-energy data set extends the kinematic region in W 1P into the transition region 



between previous diffractive J/ip measurements at HERA and fixed target experiments [2627|. 
The elastic and proton-dissociative cross sections as functions of t and W 1P are subjected to 
phenomenological fits, together with previous HI data |4||24j, and are compared with QCD 
based dipole models fl4j. 



2 Experimental Method 

2.1 Kinematics 

The kinematics of the processes ep — > e J/ip X, where X = p or Y (depicted in figure [T]), are 
described by the following variables: the square of the ep centre-of-mass energy s = (P + k) 2 , 
the square of the ^p centre-of-mass energy W 2 p = (q + k) 2 , the absolute value of the four- 
momentum transfer squared at the lepton vertex Q 2 = —q 2 — — (k — k') 2 and of the four- 
momentum transfer squared at the proton vertex t = (P — P') 2 . The four-momenta k, k', P, P' 
and q refer to the incident and scattered beam positron, the incoming and outgoing proton (or 
dissociated system Y) and the exchanged photon, respectively. 

In the limit of photoproduction, i.e Q 2 — > 0, the beam positron is scattered at small angles 
and escapes detection. In this regime the square of the 7p centre-of-mass energy can be recon- 
structed via the variable W 2 ec = s y rec , where y rec is the reconstructed inelasticity, measured 
as y rec = (Ej/tp —p z ,j/i>) / (2E e ). Here, Ej/^ and p z ,j/i> denote the reconstructed energy 
and the momentum along the proton beam direction (z-axis) of the J/ip meson and E e is the 
positron beam energy. Furthermore, the variable t can be estimated from the transverse momen- 
tum of the J/ip in the laboratory frame via the observable t rec = —p%ji^- The reconstructed 
variables W lp ^ rec and t rec are only approximately equal to the variables W 1P and t, due to their 
definition and due to the smearing effects of the detector. In particular, — p 2 T j,< is systemati- 
cally larger than t for events with a value of Q 2 close to the upper boundary of 2.5 GeV 2 used 
in the analysis. In such events the J/ip recoils against the scattered beam positron in addition 
to the proton. The measurement presented here corrects for this recoil effect by the unfolding 
procedure described below. 



2.2 Monte Carlo models 

Monte Carlo (MC) simulations are used to calculate acceptances and efficiencies for triggering, 
track reconstruction, event selection, lepton identification and background simulation. The elas- 
tic and proton-dissociative J/ip signal events are generated using the program DIFFVM [J28|, 
which is based on Regge theory and the Vector Dominance Model p9| . For J/ip production 
with proton dissociation a mass dependence of da/dM Y oc J '(My) My is implemented in 
DIFFVM. Here f{M Y ) = 1 for M Y > 3.6 GeV 2 , whereas for lower values of M Y the produc- 
tion of excited nucleon states is taken into account explicitly. The description of the forward 
energy flow and the simulated W 1P and t dependences are improved by weighting the MC sam- 
ples in W JP , t and M Y according to a functional behaviour motivated by the triple pomeron 



model [30 1 for the proton-dissociative case. The reweighting model contains seven parame- 
ters, which are adjusted to the data pTJ . QED radiation effects, which are particularly relevant 
for J/ip — > e + e~ decays, are simulated with the program PHOTOS p2| . The non-resonant di- 
lepton background is estimated using the GRAPE generator p3| , which simulates electroweak 
processes ep — > eX £ + £~ . Possible interference effects between di-lepton production via elec- 
troweak processes and J/ip decays are ignored. 



For all MC samples detector effects are simulated in detail with the GEANT program [34|. 
The MC description of the detector response, including trigger efficiencies, is adjusted using 
comparisons with independent data. Beam-induced backgrounds are taken into account by 
overlaying the simulated event samples with randomly triggered events. The simulated MC 
events are passed through the same reconstruction and analysis software as is used for the data. 



2.3 Detector 



The HI detector is described in detail elsewhere [35 36 1. Only those components essential 
for this analysis are described here. The origin of the right-handed HI coordinate system is 
the nominal ep interaction point, with the direction of the proton beam defining the positive 
z axis (forward direction). Transverse momenta are measured in the x-y plane. Polar ($) and 
azimuthal ((f)) angles are measured with respect to this frame of reference. 

In the central region (15° <i3< 165°) the interaction point is surrounded by the central 
tracking detector (CTD). The CTD comprises two large cylindrical jet chambers (CJC1 and 
CJC2) and a silicon vertex detector p7) . The CJCs are separated by a further drift cham- 
ber which improves the z coordinate reconstruction. The CTD detectors are arranged con- 
centrically around the interaction region in a uniform solenoidal magnetic field of 1.16 T. 
The trajectories of charged particles are measured with a transverse momentum resolution of 
&{Pt)/pt ~ 0.2%p T /GeV © 1.5%. The CJCs also provide a measurement of the specific 
ionisation energy loss dE/dx of charged particles with a relative resolution of 6.5% for long 
tracks. 



The liquid argon (LAr) sampling calorimeter [38 1 surrounds the tracking chambers and has 
a polar angle coverage of 4° < $ < 154°. It consists of an inner electromagnetic section with 
lead absorbers and an outer hadronic section with steel absorbers. Energies of electromag- 



netic showers are measured with a precision of a(E)/E = 12%/y/E/ GeV © 1% and those of 



hadronic showers with a(E)/E = 50%/\/E/ GeV © 2%, as determined in test beam experi- 



ments |39 40 1. In the backward region (153° <$< 178°), particle energies are measured by a 



lead-scintillating fibre spaghetti calorimeter (SpaCal) [36]. 

The calorimeters are surrounded by the muon system. The central muon detector (CMD) is 
integrated in the iron return yoke for the magnetic field and consists of 64 modules, which are 
grouped into the forward endcap, the forward and backward barrel and the backward endcap 
and cover the range 4° < •§ < 171°. 

Two sub-detectors situated in the forward direction are used in this analysis. These are 
the PLUG calorimeter, which is situated at z = 4.9 m, and consists of four double layers of 
scintillator and lead absorber, and the z = 28 m station of the forward tagging system (FTS), 
which comprises scintillator counters situated around the beam-pipe. 

HI has a four- level trigger system. The first level trigger (LI) is based on fast signals from 
selected sub-detector components, which are combined and refined at the second level (L2). The 
third level (L3) is a software based trigger using combined LI and L2 trigger information. After 
reading out the full event information events are reconstructed and subjected to an additional 
selection at a software filter farm (L4). The data used for this measurement were recorded 
using the Fast Track Trigger (FTT) [HTJ which, based on hit information provided by the CJCs, 
reconstructs tracks with subsequently refined granularity at the first two trigger levels, first in 
the x-y plane at LI and then in three dimensions at L2. 

For the data set taken at y/s ps 318GeV the luminosity is determined from the rate of 
the elastic QED Compton process ep — >■ cyp, with the positron and the photon detected in 
the SpaCal calorimeter, and the rate of Deep-Inelastic Scattering (DIS) events measured in the 
SpaCal calorimeter [42] . For the data set taken at \/s ~ 225 GeV the luminosity determination 
is based on the measurement of the Bethe-Heitler process ep — > ejp where the photon is 
detected in a calorimeter located at z = —104 m downstream of the interaction region in the 
electron beam direction. 



2.4 Event selection 

The measurement is based on two data sets, both recorded with a positron beam energy of E e = 
27.6 GeV. The first data set was taken in the years 2006 and 2007, when HERA was operated 
with a proton beam energy of 920 GeV, resulting in a centre-of-mass energy of y/s ^318 GeV. 
It corresponds to an integrated luminosity of C = 130pb _1 . The second data set was recorded 
in the last months before the shutdown of HERA in 2007, when the proton beam had a reduced 
energy of 460 GeV, resulting in \/s ~ 225 GeV. This data set corresponds to an integrated 
luminosity of C = 10.8 pb _1 . These two samples will be referred to as high-energy (HE) and 
low-energy (LE) data sets in the following. 

Photoproduction events are selected by requiring the absence of a high energy electromag- 
netic cluster, consistent with a signal from a scattered beam positron in the calorimeters. Events 
with positrons detected in the SpaCal or LAr calorimeter with with energy above 8 GeV are 
rejected. This limits the photon virtuality to Q 2 < 2.5 GeV 2 , resulting in a mean virtuality of 
(Q 2 ) = 0.1 GeV 2 . 
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The triggering of events relies on the online reconstruction of exactly two oppositely charged 
tracks with transverse momenta p T > 0.8 GeV by the FTT. This condition is verified offline 
using reconstructed tracks based on the full CTD information in the polar range 20° < $ < 
165°. 



Electrons from J/tp decays are identified using an electron estimator D [43 1, which is based 
on energy deposits and shower shape variables in the LAr calorimeter and the specific ionisation 
energy loss dE/dx measured in the CJCs. The estimator is defined such that D = 1 for genuine 
electrons and D = for background from pions. The selection of J/ip — > e + e~ events is 
performed by requiring a well identified electron with D > 0.8 in the polar range 20° < $ e < 
140°, and by observing a specific ionisation loss of the second track compatible with the electron 



hypothesis [31 1. 



In the selection of J/ip — > ^ + [i~ events one muon candidate is identified either in the 
calorimeter or in the muon system in the polar angle range of 20° < i?^ < 162.5° OTJ. In 
order to reject misidentified J/ip — > e + e~ events in this sample, the measured dE/dx values 



of both tracks must be incompatible with the electron hypothesis [31 1. The signature of a 
J/ip — > /U + /i~ event can also be mimicked by a muon from a cosmic shower passing through 
the detector. The corresponding background is rejected by an acollinearity cut and a cut on the 
timing information from the CTD pTJ . 

In order to suppress remaining non-ep background, the event vertex, which is reconstructed 
from the charged tracks in the event, is required to be within 35 cm of the nominal interaction 
point. 

The summed squared energies of the SpaCal and LAr calorimeter clusters not related to 
the J/ip decay and above 400 MeV have to satisfy the condition J2 { Ef < 2.5 GeV 2 . This 
requirement reduces the remaining background from proton-dissociative J/ip production with 
My > 10 GeV to less than 2% and from inelastic J/ip production to the per-mille level pTJ. 

The di-lepton invariant mass distributions as reconstructed from the tracks for the muon 
and the electron selection are shown in figure [2] for both the HE and LE samples. In all dis- 
tributions the J/ip peak at m« w 3.1 GeV is clearly visible. The prominent tail of the mass 
peak in the J/ip — > e + e~ channel towards low values of m ee is due to QED radiation losses 
and bremsstrahlung from the electrons, reducing their momenta. There is also background 
from non-resonant QED processes ep — > eX £ + £~ . Non-resonant diffraction contributes as a 
background to the muon channel due to pions misidentified as muons. In contrast, the elec- 
tron channel has negligible pion contamination near the J/ip mass peak due to the superior 
background rejection of the electron selection. 



2.5 Signal determination 

2.5.1 J/ip — )► ti+fj,- 

For the muon decay channel the number of reconstructed J/ip mesons is obtained from the 
invariant di-muon mass distributions m^ in bins of t rec and W 1PyTec . This is done by fitting the 



sum of a Student's t-function describing the signal and an exponential distribution for the non- 



resonant background with an extended binned log-likelihood fit using the RooFit package [44]. 
The fit model has the form 

f(N sig ,N BG ,m^;n,a,n,c) = N sig p sig (m^; /i,cr,n) + N BG p BG (m^; c) (1) 

with free shape parameters fx, a,n and c describing the probability density functions of the J/ip 
signal p sig and of the background p BG . The number of signal and background events are given 
by iVsig and N BG , respectively. The probability density functions are defined as 

/ r 2\"0.5(n+l) 

Psig(«V/J = ™ s i g f H J , r = (m w - [i)l a, n > and (2) 

Pbg("V/J =n BG e- c - m ^. (3) 

The factors n sig and n BG are chosen such that the probability densities are normalised to one 
for both p = p sig and p = p BG , in the fit range 2.3 GeV < m w < 5GeV. The small ip(2S) 
contribution is also included in the fit, modelled by a Gaussian. 

The results of the fits to the di-muon samples are shown together with the data in fig- 
ure g The fit yields to 29931 ± 217 J/ip -> //+ fi~ events for the HE data set and 2266 ± 56 
J/ip — > yU + /i _ events for the LE data sets. 



2.5.2 J/-0 -► e+e" 

For the electron decay channel the signal is determined from the invariant di-electron mass 
distributions obtained in bins of t rec and W lp ^ rec . To reconstruct the number of J/ip mesons, 
a different procedure from that used in the muon channel is employed, which minimises the 
sensitivity to details of the large radiative tail of the J/ip mass peak visible in figure [2] and dis- 
cussed above. The non-resonant background, modelled using the QED process ep — ¥ eX £ + £~ 
as simulated with GRAPE, is subtracted from the data. This is possible due to the negligi- 
ble contamination from particles other than electrons at and above the J/ip mass peak. The 
normalisation of the simulated QED background is determined prior to the background sub- 
traction by fitting the background to the overall invariant mass distribution in the mass window 
3.75 < m ee < 5 GeV above the ip(2S) mass, where only the QED contribution is present. 
Within errors this normalisation factor is consistent with unity. 

After background subtraction the remaining events are counted within a window of 2.3 < 
m ee < 3.3 GeV around the nominal J/ip mass peak. This yields 23662 ± 177 J/ip — ¥ e + e~ 
events for the HE data set and 1760 ± 47 J/ip — > e + e~ events for the LE data sets. These 
numbers of events are then corrected to account for the fraction of signal events outside the 
counting window, which is close to 5% as determined using the J/ip MC simulation. Within 
the counting window the J/ip MC simulation describes the behaviour of the radiative tail in the 
data well. 



2.6 Experimental signatures of elastic and proton-dissociative processes 

Proton-dissociative candidate events are identified by requiring either a large value of \t rec \ > 
1.5 GeV 2 or energy deposits in the HI forward detectors, caused by fragments of the proton- 
dissociative system. Three subdetectors, situated at different locations, are used in this analysis 
to measure activity in the forward direction, using the following requirements. 

• At least one cluster well above the noise level is found in the forward part of the LAr, 
with an energy above 400 MeV and ■§ < 10°. 

• The summed energy of all clusters in the PLUG calorimeter is above 4 GeV, where all 
clusters above the threshold level of 1.2 GeV are considered. 

• Activity is observed in at least one scintillator of the FTS station situated at z = 28 m. 

If at least one of these conditions is fulfilled, the event is flagged as tagged. Identical tagging 
methods are applied in the e + e~ and /i + fi~ channels. 

In figure[3]the simulated tagging efficiencies and tagging fractions observed in data and sim- 
ulation are shown as functions of W JPtrec and — t rec . The tagging fractions are obtained from the 
e + e sample, and contain QED contributions in addition to di-electron events from diffractive 
J/ip production. In order to enrich it with genuine J ftp decays, the sample is restricted to in- 
variant masses in the window m ee = 2.3 — 3.3 GeV. The tagging fractions observed in the data 
are compared to the simulation. The simulation is based on the MC generators DIFFVM for 
elastic and proton-dissociative J/ip production and GRAPE, which is used to describe the QED 
background. The uncertainty in the simulation due to the tagging of the forward energy flow is 
represented by the shaded bands. The tagging efficiency and fraction show a flat behaviour as a 
function of W 7Ptrec . A steep rise of the tagging fraction is observed as a function of t rec , which 
reflects the relative elastic and proton-dissociative contribution in data. 

An unambiguous event-by-event distinction between elastic and proton-dissociative events 
is not possible with the HI detector. Proton-dissociative events can be misidentified as elastic 
events if the outgoing dissociated proton remains undetected due to the limited acceptance of 
the forward detectors. On the other hand, elastic events may have significant energy deposits 
in the forward detectors due to possible beam induced background and may be misidentified 
as proton-dissociative events. However, since the forward energy flow is modelled by the MC 
simulation, elastic and proton-dissociative cross sections can be unfolded on a statistical basis. 



2.7 Unfolding 

Regularised unfolding is used to determine the elastic and proton-dissociative cross sections 
in bins of t and W 1P from the number of events observed as a function of t rec and W 7Pirec , 
respectively, and from the tagging information as described in the previous subsection. The 
general procedure is described in [43} 45 46 1 and the references therein. In the following only 



the aspects most relevant to this analysis are summarised; further details are discussed in pT| . 
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All efficiency corrections and migration effects are described by a response matrix A, which 
correlates the number of reconstructed J/ip events in each analysis bin, represented by the 
vector y TCC , with the true distribution x tTUC via the matrix equation y TCC = AaJtrue- The matrix 
element Ay gives the probability for an event originating from bin j of x truc to be measured in 
bin i of y rec . The unfolded "true" distribution is obtained from the measured one by minimising 
a x 2 -function x 2 (#tru e ; 2/ rec ) by variation of x tIUC , with a smoothness constraint determined by 
a regularisation parameter. This parameter is chosen such that the correlations in the covariance 
matrix of the unfolded distribution x tmc are minimised. 

Two types of response matrix A are used: one to unfold differential cross sections as a 
function of t, and one to unfold differential cross sections as a function of W 1P . The response 
matrices are calculated from the simulation and are defined such that the elastic and proton- 
dissociative differential cross sections are determined simultaneously. By using the tagging 
information for small values of \t rec \ < 1.5 GeV 2 , the elastic and proton-dissociative cross 
sections are disentangled. Since the region of large values of \t rec \ is completely dominated 
by proton dissociation, no tagging condition is applied. Further, two reconstructed bins are 
associated with each bin at the truth level, in order to provide sufficiently detailed information 
on the probability distribution and to improve the accuracy of the unfolding procedure. 

The unfolding procedure is applied separately for the HE and the LE data sets. The response 
matrices for the LE data set are similar to those for the HE case. However they contain fewer 
bins due to the smaller number of events. 

In figure [4] control distributions are shown for W lp ^ rec and — t rec separately for the fi + fi~ 
sample and the e + e~ sample. Both samples are restricted in meg to the J/ift peak region, which 
is chosen for the /U + /i~ sample to be 2.8 < m w < 3.3 GeV. For the e + e~ sample this region 
is enlarged to 2.3 < m ee < 3.3 GeV in order not to cut into the radiative tail. The relative frac- 
tions of the elastic and proton-dissociative events simulated with DIFFVM as determined in the 
unfolding procedure, are also shown in figure |4j The contribution from the ip(2S) resonance 
is taken from the simulation, normalised using a previous measurement [ |20[ . For the e + e~ 
sample, the QED background simulated with GRAPE is indicated and normalised as described 
above. For the control distributions of the /i + fi~ sample the non-resonant background is sub- 
tracted from the data using a side band method OTJ. This background contains a contribution 
from non-resonant diffractive events, due to pions misidentified as muons, in addition to the 
QED background. The data in all distributions are well described by the simulation. 

2.8 Cross section determination and systematic uncertainties 

The cross sections are measured for the kinematic ranges as defined in table [TJ From the un- 
folded number of events in each signal bin i for the reaction 'jp — > J/ip — > ££, the bin-averaged 
cross sections are obtained as 

dt " <&t C ■ B{££) ■ Mi ' 

and 



-1 
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Data Set 


E p 


Process 


Q 2 


My 




\t\ 


w lp 


HE 


920 GeV 


elas 
pdis 


< 2.5 GeV 2 


m p 
m p - 10 GeV 




< 8 GeV 2 


40 -110 GeV 


LE 


460 GeV 


elas 
pdis 


< 2.5 GeV 2 


m p 
m p - 10 GeV 


< 


5«, 8 GeV 2 


25 - 80 GeV 



Table 1: Kinematic range of the analysis. The phase space for elastic and proton-dissociative 
J/ip processes is indicated by elas and pdis, respectively. The high- and low-energy data sets 
are denoted by HE and LE. W The phase space restriction is applied only for the da/dt cross 
section measurement. 



where the variable $^ is the transverse polarised photon flux [28 1, $ 7 A 7P the transverse po- 
larised photon flux per W 1P bin, Atj the bin width in t, 11 = ee or /x/x depending on the decay 
channel, N i>t> u and N i>w « are the numbers of unfolded signal events in the corresponding 
bins of t or W 1P , C is the integrated luminosity, and 13 (ee) = 5.94%, £>(/x/x) = 5.93% are the 



J/ip branching fractions [47|. 



The systematic uncertainties on the J/ip cross section measurement are determined by im- 
plementing shifts due to each source of uncertainty in the simulation and propagating the re- 
sulting variations in the unfolding matrices to the result. Those uncertainties which are uncorre- 
cted between the two decay modes are classified as individual systematic uncertainties, while 
the uncertainties correlated between the e + e~ and ii + /x~ samples are referred to as common 
systematic uncertainties. 

The individual systematic uncertainties are as follows. 

Lepton identification The efficiency of the simulated muon identification is reweighed to agree 
with that determined from data. The efficiency was determined with a J/ip — ¥ /x + /x~ 
sample, selected with at least one identified muon. The second muon is then probed to 
evaluate the single muon identification efficiency. The uncertainty on these weights is 



determined from the remaining difference between the simulation compared to data [31 1. 
The resulting uncertainty on the cross sections is 2% at most. 

The cut value on the electron discriminator D is varied by ±0.04 around its nominal value 
of 0.8, which covers the differences in the /J-distribution between simulation and data. 
The uncertainty propagated to the cross section is below 2%. 

Signal extraction The uncertainty on the number of signal events due to the fitting proce- 
dure of the m^ invariant mass distributions is determined by a bias study as described 



in [31 1 and is typically ss 1% but can rise to m 5% for the lowest W 1P bin of the proton 



dissociative cross section. 

The uncertainty on the background subtraction procedure for the e + e~ sample is esti- 
mated by determining the background normalisation factor with data at very low invariant 
di-electron masses m ee and agreement with the default method is found within 20%. The 
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corresponding variation on the background is propagated to the differential cross sections 
which vary between 3% for bins with a low background to 11% for bins with a higher 
background contribution. 

Branching ratio The relative uncertainty on the branching ratio for the muon and electron 
decay channels is 1% [ |47| . 

The following systematic uncertainties are have components contributing to the channel- 
specific individual and the common systematic uncertainties. 

Trigger The trigger efficiency is typically 80% and is taken from the simulation. The trig- 
ger simulation is verified by a comparison to data in a sample of J/ip mesons in deep- 
inelastic-scattering triggered independently on the basis of the scattered beam positron. 
A small difference of 3% is observed between the data and the simulation for J/ift events 
decaying into muons. This difference is accounted for by a corresponding upwards shift 
of the efficiency in the simulation. No such correction is necessary for electrons. The 
remaining uncertainty is estimated to be 2% uncorrected between the e + e~ and fi + ji~ 
samples, i.e. treated as individual uncertainties, and 2% correlated between the two de- 
cay channels, i.e. treated as a common uncertainty. 

Track finding efficiency The uncertainty due to the track reconstruction efficiency in the CTD 
is estimated to be 1% per track [ |48| . For electron tracks an additional 1% is applied, 
to account for the different hit finding efficiency due to bremsstrahlung effects. Since 
the uncertainty on the track finding efficiency affects both selected tracks coherently, a 
common uncertainty of 2% is applied to both samples and an additional 2% is applied for 
the electron sample. 

The following common systematic uncertainties are considered. 

Tagging The systematic uncertainty arising from the tagging condition is estimated by varying 
separately the simulated tagging efficiency for each detector used. The variations cover 
any possible shift in the individual relative efficiency distributions, and are 20% for the 



condition from the forward LAr calorimeter, 5% for the PLUG and 1% for the FTS pi] . 
The resulting uncertainties on the cross sections are typically a few percent, but reach 
30% at the highest \t\ values of the elastic da/dt cross section. 

Empty calorimeter The uncertainty on the cut ensuring an empty calorimeter is obtained by 
varying the maximum allowed ^^f fr° m 2.25 GeV 2 to 2.75 GeV 2 in the simulation. 
This results in an uncertainty of typically 5% for the proton-dissociative cross sections. 
For the elastic cross sections this variation is negligible for most bins, except for the 
highest bin in \t\, where it reaches to 13%. 

MC modelling The model uncertainty in the MC simulation due to uncertainties in the depen- 
dences on t, W 1P and M Y is determined by varying the fit parameters of the weighting 



procedure within the errors obtained in a dedicated fit of the forward energy flow [ 3 1 1 
For the cross section as a function of W 1P the corresponding uncertainties are below 



'/o, 



whereas for the cross sections differential in t, values around 10% are obtained for the 
high \t\ bins. 
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Luminosity The integrated luminosity is known to within ±2.7% for the HE data set and to 



within ±4% for the LE data set [42 1. 



ip(2S) background Background from ip(2S) decays to J/ip X is estimated to contribute 4% 
to the selected J/ip events, and is subtracted from the data prior to the unfolding proce- 
dure p0|. The cross section measurements are affected by an uncertainty of 1.5%. 



Q 2 dependence The Q 2 dependence of the cross section is parametrised as a rp oc (m\ + Q 2 ) 
Q. The corresponding systematic uncertainty is obtained by varying the parameter n in 
the range 2.50 ± 0.09. The cross sections are affected by less than 1%. 

The differential cross sections obtained from the electron and the muon data agree within 
uncertainties. The two measurements are combined by taking into account their individual 
uncertainties. This combination procedure involves the numerical minimisation of a standard 
X 2 function including the full statistical error matrix and the correlated systematic errors with 



nuisance parameters, similar to that defined in [49 50 1. All individual uncertainties are in- 
corporated within this procedure, whereas the common uncertainties are considered after the 
combination only. The consistency of the data sets can be verified by looking at the result- 
ing nuisance parameters. None of the nuisance parameters shifts by more than one standard 
deviation. 

Figure [5] shows the result of the combination for the elastic and proton-dissociative cross 
sections as a function of W yp . The input data obtained in the electron and muon decay channels 
are shown together with the combined data. 



3 Results 

The elastic and proton-dissociative differential J/ip cross sections as functions of t and W 7P are 
measured in the kinematic ranges defined in table[I]using the decay channels J/ip — > ^i + ^T and 

J/ip ->■ e + e~. 

Tables |4}|5]and|6]list the combined data points for all cross sections together with their uncer- 
tainties and all common systematic uncertainties. The input data to the combination procedure, 
including all individual systematic uncertainties together with the full covariance matrices of 



the combined results can be found in [51 1. 



3.1 t dependence 

Figure [6] shows the measured elastic and proton-dissociative cross sections differential in —t, 
separately for the LE and HE data sets. The cross sections fall steeply with increasing —t, and 
shows a clear difference between the shapes of the proton-dissociative and elastic distributions. 
The proton-dissociative cross section levels off for very low values of \t\ . There is a phase space 
effect such that for small \t\ it is not possible to produce large masses of M Y . 

In figure [7] the proton-dissociative measurement from the HE data set as a function of —t 
is compared to a previous analysis p4~} covering the region of high |t|, which is completely 
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dominated by proton-dissociative events. The high \t\ data [24| are adjusted to the W 1P , Q and 
My ranges of the present analysis by applying a phase space correction of about 7%. Comparing 
the two measurements, the present proton-dissociative cross sections extend the reach to small 
values of \t\. In the overlap region 2 < \t\ < 8 GeV 2 the two measurements agree. 

The elastic and proton-dissociative differential cross sections da/dt are fitted simultane- 
ously, using a x 2 -function [ |49lp0| based on the error matrix obtained in the combination pro- 
cedure and all common systematic uncertainties. The elastic cross section is parametrised as 
da/dt = N e i e~ bel ^. For the proton-dissociative cross section da/dt = N pd (1 + (6 pd /n)|£|)~" 
is chosen, which interpolates between an exponential at low \t\ and a power law behaviour at 
high values of \t\. The fits are performed separately for the HE and the LE measurements. In 
the case of the HE data the previously measured high \t\ data are included in the fit. This fit 
yields a value of x 2 /NDF = 26.6/18 after excluding the two lowest t data points in both the 
elastic and the proton-dissociative channel. For fit of the LE data set, the parameter n is fixed to 
the value obtained from the HE data set, since the LE data are not precise enough to constrain 
b p d and n simultaneously. The obtained parametrisations for the elastic and proton-dissociative 
cross sections are compared to the data in figure [6] and figure [7] Table [2] summarises the fit 
parameters and their uncertainties. 

The elastic cross section data for —t > 0.1 GeV are well described by the exponential 
parametrisation. They fall much faster with increasing \t\ than the proton-dissociative cross 
section even at small \t\, which is reflected in the values for b e i and b p d- The value extracted 
for b e i is compatible with previous results Q, although the previous fit was done as a function 
of p 2 T j,, rather than —t. Some difference between the b el values for the LE and HE data is 
expected [4] due to the different ranges in W 1P corresponding to (W 1P ) = 78 GeV for the HE 
data set and (W 1P ) = 55 GeV for the LE data. 



3.2 Energy dependence 

The measured elastic and proton-dissociative cross sections as a function of W 7P are shown in 
figures [8} The elastic and proton-dissociative cross sections are of similar size at the lowest 
W JP = 30 GeV accessed in this analysis. The elastic cross section rises faster with increasing 
W 1P than the proton-dissociative one. The ratio of the proton-dissociative to the elastic cross 
section as a function of W 1P is also shown in figure [8j The ratio decreases from 1 to 0.8 
as W JP increases from 30 GeV to 100 GeV. When calculating the ratio no attempt is made to 
extrapolate the elastic measurement to — t = 8 GeV 2 . The corresponding correction is estimated 
to be smaller than 1%. 

In figure [9] the elastic cross section measurements of this analysis are compared to previ- 
ous measurements at HERA [[4l[5|. The LE data extend the range accessible in W rp to lower 
values when compared to previous HI measurements Q. The HE data have a large overlap 
with previous HI measurements in the region 40 GeV < W lp < 110 GeV and show a similar 
precision. Within normalisation uncertainties, the previous measurements and the new data are 
in agreement. 
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Data period Process Parameter Fit value Correlation 



HE -fp^J/ipp b d (4.88 ± 0.15) GeV 



-2 



p(b e i,N el ) = 0.50 
p(bel,b pd ) =0.49 
p(b el ,n) =-0.21 
p(b e i,N pd )=0.6S 



p(N e i,b pd ) =0.23 

N e t (305 ± 17) nb/GeV 2 rt^d.") =-°- 07 

p(N el ,N pd )=0M 



p(6 pd ,n) =-0.78 
p(b pd , N pd )=0.16 



1P -> J/V F 6pd (1.79 ± 0.12) GeV" 2 

n 3.58 ±0.15 P (n, Jv pd )=-o.46 

iVpd (87 ± 10) nb/GeV 2 

/9(6 e! ,Ar ei ) = 0.37 

LE ^p^J/ipp b d (4.3±0.2)GeV~ 2 P (b el ,b pd ) =o.w 

p(b el ,N pd )=0Al 

N el (213 ± 18) nb/GeV 2 $££££? 

7P -> J/V y 6 pd (1.6 ± 0.2) GeV" 2 p(fe pd , w pd )=o.53 

n 3.58 (fixed value) 

A^ d (62 ± 12) nb/GeV 2 

Table 2: Parameter values obtained from the fits to the differential cross sections da/dt, includ- 
ing their errors and correlations. The fit functions are described in the text. HE and LE denote 
the high- and low-energy data sets, respectively. 

The measured elastic and proton-dissociative cross sections as a function of W 1V , shown 
in figure [8j are fitted simultaneously, taking into account the correlations between the proton- 
dissociative and the elastic cross sections. The fit also includes data from a previous measure- 
ment [4 1 shown in figure [8} with a normalisation uncertainty of 5% and all other systematic 
uncertainties treated as uncorrelated. As parametrisation two power law functions of the form 
a = N (W lp /W 1P fl) with W 1P $ = 90 GeV are used with separate sets of parameters for the 
elastic and the proton-dissociative cases. The % 2 -function is defined in the same manner as for 
fits of the t-dependences. 

The result of the fit is compared to the measurements in figures [8] and in figures |9j The 
parametrisation describes the data well (x 2 /NDF = 32.6/36). The fitted parameters are given 
in table [3] together with their uncertainties and correlations. In Regge phenomenology the pa- 
rameter 5 can be related to the pomeron trajectory a(t) = a(0) + a' • t by S(t) = 4(a(t) — 1). 
Using the values a' el = 0.164 ± 0.028 ± 0.030 GeV" 2 El and a' pd = -0.0135 ± 0.0074 ± 
0.0051 GeV -2 J24}, together with the mean values of t for the elastic and proton-dissociative 
measurements, (t) = —0.2 GeV 2 and (t) = —1.1 GeV 2 , one can estimate a(0) for the elas- 
tic and proton-dissociative process from these measured parameters. The obtained values of 



a(0) e i = 1.20 ± 0.01 and a(0) pd = 1.09 ± 0.02 are in agreement with the results from [4,5 52]. 

The direct comparison between 5 e i and 5 p d is made by looking at the ratio of the two cross 
sections, shown in figure p3J The ratio is parametrised as N R (W 7P /W 7Pi o) R with W 7P)0 = 
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90GeV, N R = N pd /N el = 0.81 ± 0.10 and S R 



> P d 



-0.25 ± 0.06, taking all 



correlations into account. Qualitatively the decrease of this ratio with increasing W 1P has been 
predicted in [[8| as a consequence of the non-unit and W JP dependant survival probability for 
the proton dissociation process. 



In figure 10 a compilation of cross section measurements for the elastic J/ip cross section 
is shown as a function of W™. The LE data from the present analysis close the gap to data 
from fixed target experimentsH[26 27 1 at low W 1P . The fixed target data exhibit a lower nor- 



malisation and a steeper slope than observed at HERA. Also shown are recent results from the 



LHCb experiment [53 1. The extrapolated fit function for the elastic J/ip cross section is able to 
describe the LCHb data points at high W 1P well. 



W, 



Following |14j the obtained value of 5 can for large photon-proton centre-of-mass energies, 
rp ^$> mj/ip, be related to a leading-order gluon-density parametrised as x-g(x,fj 2 ) = N-x~ x 
via 5 e i ~ 4 ■ A. The scale of J/ip photoproduction is often taken to be /i 2 = 2.4 GeV 2 . The 
observed value \j/$ = 0.168 ± 0.008 is in remarkable agreement with \i n d(Q 2 = 2.5 GeV 2 ) = 
0.166 ± 0.006 obtained from fits to inclusive DIS cross sections [49] . Skewing effects [ 14[54| 
are ignored in this comparison. 



In [ 14 1 both a leading order and a next-to-leading order gluon-density are derived, via fits to 
previous J/ip measurements at HERA |4} [2l] - |23| . The fit results obtained in [ 14 1 are compared 
with the data in figure 1 1 Both fits are also extrapolated from the W 1P range of the input data 
to higher W 1P and compared with the LHCb measurement. The leading-order fit describes the 
LHCb data well, whereas the next-to-leading order fit lies above the LHCb cross sections. 



Process 


Parameter 




Fit value 


Correlation 


7P — >■ J/ipp 


6 e , 






0.67 ±0.03 


p(S e i,N el )=-om 

p(Sel,S pd ) =0.01 

p(5 e i,N pd )=0M 




N el 






81 ± 3 nb 


p(N el ,S pd ) =-0.27 
p(N el ,N pd )=-0.U 


7P — > J/ip Y 


5 P d 
N pd 






0.42 ± 0.05 
66±7nb 


p(5 pd ,N pd )=0.09 


Ratio 


5r = 


8pd - 


- 5 e l 


-0.25 ±0.06 


p(S r ,N R )=0.l4 




N R ~- 


= N pd 


/N el 


0.81 ±0.11 





Table 3: Parameter values obtained from the fit to the cross sections as a function of W^ n , in 



IV 



eluding their errors and correlations. The fit functions are described in the text. The parameters 
for the ratio of the two functions are also given. 



'The data from |26| and |27| have been updated using recent measurements of branching ratios (47). The data 
from [26 1 are also corrected for contributions from inelastic processes, see |51 1 for more details. 
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4 Summary 

Photoproduction cross sections for elastic and proton-dissociative diffractive J/ip meson pro- 
duction have been measured as a function of t, the four-momentum transfer at the proton ver- 
tex, and as a function of W 7P , the photon proton centre-of-mass energy in the kinematic ranges 
|t| < 8 GeV 2 , 25 GeV < W 7P < 110 GeV and for the proton-dissociative case M Y < 10 GeV. 
The data were collected in positron-proton collisions with the HI detector at HERA, at a centre- 
of-mass energy of y/s w 318 GeV and \/s « 225 GeV. Measurements in the electron and 
muon decay channels are combined, and are parametrised using phenomenological fits. 

The elastic and the proton-dissociative cross sections are extracted simultaneously. Using 
this technique, a precise measurement of proton-dissociative J ftp production was performed in 
the range of small \t\ for the first time. The data taken at low centre-of-mass energies close the 
gap between previous HI measurements and fixed target data. 

The data agree well with previous HERA measurements and with a model based on two 
gluon exchange. The W^p-dependence of the proton-dissociative channel is found to be signifi- 
cantly weaker than that of the elastic channel. 



5 Acknowledgements 

We are grateful to the HERA machine group whose outstanding efforts have made this ex- 
periment possible. We thank the engineers and technicians for their work in constructing and 
maintaining the HI detector, our funding agencies for financial support, the DESY technical 
staff for continual assistance and the DESY directorate for support and for the hospitality which 
they extend to the non-DESY members of the collaboration. We would like to give credit to 
all partners contributing to the WLCG computing infrastructure for their support for the HI 
Collaboration. 



References 

[1] P. D. B. Collins, "An Introduction to Regge Theory and High Energy Physics," Cambridge 
University Press, 1977. 

[2] A. Donnachie and P. V. Landshoff, Phys. Lett. B 348 (1995) 213 [hep-ph/94TT36§. 

[3] A. Levy, "Low x physics at HERA," Lecture Notes in Physics, 496, (1997), 347-477, 
DESY-97-013, TAUP-2398-96. 



[4] A. Aktas et al. [HI Collaboration], Eur. Phys. J. C 46 (2006) 585 [hep-ex/05 10016 1. 

[5] S. Chekanov etal. [ZEUS Collaboration], Eur. Phys. J. C 24 (2002) 345 [hep-ex/0201043]. 

[6] V. A. Khoze, A. D. Martin and M. G. Ryskin, Phys. Lett. B 643 (2006) 93 [hep- 



ph/0609312]. 

18 



[7] E. G. S. Luna et al, Eur. Phys. J. C 59 (2009) 1 [arXiv:0807.4115| 



[8] E. Gotsman et al, Eur. Phys. J. C 52 (2007) 295 [hep-ph/0702053]. 
[9] M. G. Ryskin, Z. Phys. C 57 (1993) 89. 
[10] S. J. Brodsky et al, Phys. Rev. D 50 (1994) 3134 [hep-ph/9402283|. 



[11] D. Y. Ivanov et al, Eur. Phys. J. C 34 (2004) 297 [hep-ph/0 401131] . 
[12] N. Brambilla etal, Eur. Phys. J. C 71 (2011) 1534 [arX iv:1010.5"827 | 
[13] T. Teubner, hep-ph/99 10329. 



[14] A. D. Martin et al, Phys. Lett. B 662 (2008) 252 [arXiv: 0709.4406]. 

[15] M. Diehl, Phys. Rept. 388 (2003) 41 [hep-ph/0307382j. 

[16] L. Frankfurt, M. McDermott and M. Strikman, JHEP 0103 (2001) 045 [hep-ph/0009086]. 

[17] A. Caldwell and H. Kowalski, Phys. Rev. C 81 (2010) 025203. 

[18] S. Aid etal [HI Collaboration], Nucl. Phys. B 472 (1996) 3 fhejvex/9603005]. 

[19] C. Adloff et al [HI Collaboration], Phys. Lett. B 483 (2000) 23 [hep-ex/0003020|. 



[20] C. Adloff etal [HI Collaboration], Phys. Lett. B 541 (2002) 251 [hep-ex/0205107]. 

[21] J. Breitweg etal [ZEUS Collaboration], Z. Phys. C 75 (1997) 215 [ |hep-ex/97040T3| |. 

[22] J. Breitweg etal [ZEUS Collaboration], Eur. Phys. J. C 6 (1999) 603 | |hep-ex/9 8080201. 

[23] S. Chekanov etal [ZEUS Collaboration], Nucl. Phys. B 695 (2004) 3 | |hep-ex/0 4040081. 



[24] A. Aktas etal [HI Collaboration], Phys. Lett. B 568 (2003) 205 [hep-ex/0306013] 



[25] S. Chekanov etal [ZEUS Collaboration], JHEP 1005 (2010) 085 [a rXiv:0910.1235| |. 

[26] M. E. Binkley et al, Phys. Rev. Lett. 48 (1982) 73. 

[27] B. H. Denby et al, Phys. Rev. Lett. 52 (1984) 795. 

[28] B. List and A. Mastroberardino, "DIFFVM - A Monte Carlo generator for diffractive 
processes in ep scattering,", Proc. of the Workshop on Monte Carlo Generators for HERA 
Physics, eds. A.T. Doyle et al., DESY-PROC- 1999-02 (1999) 396. 

[29] J. J. Sakurai, Annals Phys. 11 (1960) 1; 

M. Gell-Mann and F. Zachariasen, Phys. Rev. 124 (1961) 953; 

T. H. Bauer et al, Rev. Mod. Phys. 50 (1978) 261 [Erratum-ibid. 51 (1979) 407]. 

[30] N. P. Zotov, and V. A. Tsarev, Sov. Phys. Usp. 31 (1988) 119 



19 



[31] F. Huber, "Elastic and proton-dissociative J/ip Photoproduction at low W 1P with the HI 
Detector at HERA", Ph.D. thesis, Universitat Heidelberg (2012), DESY-THESIS-2013- 
004 (available at http://www-hl.desy.de/psfiles/theses/). 

[32] E. Barberio and Z. Was, Comput. Phys. Commun. 79 (1994) 291; 

Z. Was, P. Golonka and G. Nanava, PoS ACAT (2007) 071 [ |arXiv:0707.30 44]. 

[33] T. Abe, Comput. Phys. Commun. 136 (2001) 126 [hep-ph/00 12029]. 

[34] Geant 3, R. Brun et al, "Geant3," CERN-DD/EE/84- 1 (1987). 

[35] I. Abt et al, [HI Collaboration], Nucl. Instr. and Meth. A 386 (1997) 310; ibid, 348. 

[36] R.D. Appuhn et al, Nucl. Instrum. Meth. A 386 (1997) 397. 

[37] D. Pitzl et al, Nucl. Instrum. Meth. A 454 (2000) 334 [hep-ex/0002044]. 



[38] B. Andrieu et al. [HI Calorimeter Group], Nucl. Instrum. Meth. A 336 (1993) 460. 

[39] B. Andrieu et al [HI Calorimeter Group], Nucl. Instrum. Meth. A 350 (1994) 57. 

[40] B. Andrieu et al [HI Calorimeter Group], Nucl. Instrum. Meth. A 336 (1993) 499. 

[41] A. Baird et al, IEEE Trans. Nucl. Sci. 48 (2001) 1276 [ |hep-ex/01040 10|; 
D. Meer et al, IEEE Trans.Nucl.Sci., 49, (2002) 357; 

N. Berger et al, IEEE Nuclear Science Symposium Conference Record, volume 3, (2004) 
1976. 



[42] F. D. Aaron et al [HI Collaboration], Eur. Phys. J. C 72 (2012) 2163 [arXiv: 1205.2448] 



[43] F. D. Aaron et al [HI Collaboration], Eur. Phys. J. C 72 (2012) 2148 [arXiv: 1206.4346]. 
[44] W. Verkerke and D. P. Kirkby, eConf C 0303241 (2003) MOLT007 [jphysics/0306116]. 



[45] F D. Aaron et al [HI Collaboration], Eur. Phys. J. C 66 (2010) 17 [ arXiv : 09 10.563 lfl . 

[46] F D. Aaron et al [HI Collaboration], Eur. Phys. J. C 72 (2012) 1970 [ |arXiv:1111.05 84]. 

[47] J. Beringer et al. [Particle Data Group Collaboration], Phys. Rev. D 86 (2012) 010001. 

[48] M. Brinkmann, "Measurement of the D* ± Meson Production Cross Section and Fg at 
High Q 2 in ep Scattering at HERA", Ph.D. thesis, University Hamburg (2010), DESY- 
THESIS-2010-016 (available at http://www-hl.desy.de/psfiles/theses/). 



[49] F D. Aaron et al. [HI Collaboration], Eur. Phys. J. C 63 (2009) 625 [arXiv:0904.0929] 



[50] F D. Aaron et al. [HI and ZEUS Collaborations], JHEP 1001 (2010) 109 



[arXiv:091 1.0884]. 



[51] Additional material of this analysis is provided at URL 

http: //www-hl.desy. de/publications/Hlpublicati on.short list .html 



[52] F D. Aaron et al. [HI Collaboration], JHEP 1005 (2010) 032 [arXiv:09 10.5831] 

20 



[53] R. Aaij et al. [LHCb Collaboration], J. Phys. G 40 (2013) 045001 [arXiv: 1301.7 084] 

[54] L. Favart, M. V. T. Machado and L. Schoeffel, hep-ph/051 1069. 

[55] G. D. Lafferty and T. R. Wyatt, Nucl. Instrum. Meth., A355 (1995) 541. 



21 



W^, range (W%) 
[GeVJ LGeVJ 



[nb] [nb] 



[nb] 



L%J 



,,-Trk.o.r, 
"sys 

[%] 






4 2S <5 £ » <5 £ ' (5 LArl 

"sys u sys sys "sys 

[%] [%\ [%] [%] 



IPLUG rFTS IMC Model iQ 2 fR LT SEC 
sys sys sys sys sys sys 

[%} [%] [%] [%] [%l [%] 



High energy 


data period for elastic 


■J j'ljj production 




























40.0 - 46.5 


43.2 


0.0158 


50.7 


4.9 


2.1 


62 


2.0 


2.0 


1.5 


2.7 


-2.6 


0.6 


-0.1 


-7.1 


-0.1 


0.0 


1.4 


46.5-53.5 


50.0 


0.0144 


59.5 


5.8 


2.2 


69 


2.0 


2.0 


1.5 


2.7 


-2.6 


0.6 


-0.1 


-7.3 


-0.1 


().() 


1.4 


53.5-61.2 


57.3 


0.0131 


61.8 


6.2 


2.7 


71 


2.0 


2.0 


1.5 


2.7 


-2.6 


0.6 


-0.1 


-7.4 


-0.1 


().() 


1.4 


61.2-69.4 


65.3 


0.0120 


67.6 


6.2 


2.5 


71 


2.0 


2.0 


1.5 


2.7 


-2.6 


0.6 


-0.1 


-6.6 


-0.1 


().() 


1.4 


69.4 - 78.4 


73.9 


0.0112 


72.4 


6.4 


2.6 


71 


2.0 


2.0 


1.5 


2.7 


-2.6 


0.6 


-0.1 


-6.3 


-0.1 


0.0 


1.4 


78.4-88.0 


83.2 


0.0103 


79.9 


7.0 


3.0 


69 


2.0 


2.0 


1.5 


2.7 


-2.6 


0.6 


-0.1 


-6.0 


-0.1 


0.0 


1.4 


88.0-98.5 


93.3 


0.0096 


84.4 


7.0 


3.0 


69 


2.0 


2.0 


1.5 


2.7 


-2.6 


0.6 


-0.1 


-5.5 


-0.1 


0.0 


1.4 


98.5- 110.0 


104.3 


0.0089 


86.7 


7.3 


3.7 


65 


2.0 


2.0 


1.5 


2.7 


-2.6 


0.6 


-0.1 


-5.2 


-0.1 


0.0 


1.4 



High energy 


data period for proton 


dissociative J/4 


production 
























40.0 - 46.5 


43.2 


0.0158 


46.0 


6.0 


2.3 


54 


2.0 


2.0 


1.5 


2.7 


9.4 


-2.2 


0.5 


-3.9 


0.1 


0.0 


-4.3 


46.5-53.5 


50.0 


0.0144 


52.1 


6.5 


2.3 


61 


2.0 


2.0 


1.5 


2.7 


9.4 


-2.2 


0.5 


2.5 


0.1 


0.0 


-4.3 


53.5-61.2 


57.3 


0.0131 


58.7 


7.4 


2.3 


61 


2.0 


2.0 


1.5 


2.7 


9.4 


-2.2 


0.5 


3.5 


0.1 


0.0 


-4.3 


61.2-69.4 


65.3 


0.0120 


58.7 


7.5 


2.2 


63 


2.0 


2.0 


1.5 


2.7 


9.4 


-2.2 


0.5 


4.6 


0.1 


0.0 


-4.3 


69.4 - 78.4 


73.9 


0.0112 


61.5 


8.0 


2.4 


62 


2.0 


2.0 


1.5 


2.7 


9.4 


-2.2 


0.5 


4.8 


0.1 


0.0 


-4.3 


78.4-88.0 


83.2 


0.0103 


67.7 


8.7 


2.6 


60 


2.0 


2.0 


1.5 


2.7 


9.4 


-2.2 


0.5 


4.6 


0.1 


0.0 


-4.3 


88.0-98.5 


93.3 


0.0096 


69.8 


9.0 


2.7 


59 


2.0 


2.0 


1.5 


2.7 


9.4 


-2.2 


0.5 


4.8 


0.1 


0.0 


-4.3 


98.5-110.0 


104.2 


0.0089 


68.8 


9.0 


3.0 


54 


2.0 


2.0 


1.5 


2.7 


9.4 


-2.2 


0.5 


4.6 


0.1 


0.0 


-4.3 



Low energy data period for elastic J/ip production 



25.0 


39.0 


3 1 .9 


0.0465 


39.7 


4.9 


3.4 


62 


2.0 


2.0 


1.5 


- 4.0 


-3.4 


0.8 


-0.1 


-6.4 


-0.1 


0.0 


1.9 


39.0 


57.0 


47.9 


0.0359 


55.4 


5.6 


3.3 


64 


2.0 


2.0 


1.5 


4.0 


-3.4 


0.8 


-0.1 


-5.1 


-0.1 


0.0 


1.9 


57.0 


80.0 


68.4 


0.0284 


66.4 


6.8 


4.3 


64 


2.0 


2.0 


1.5 


4.0 


-3.4 


0.8 


-0.1 


-4.7 


-0.1 


0.0 


1.9 



Low energy data period for proton dissociative J/ip 


production 
























25.0-39.0 31.9 0.0465 42.0 8.1 


4.8 59 


2.0 


2.0 


1.5 


4.0 


12.0 


-2.9 


0.3 


-4.6 


0.1 


0.0 


-6.4 


39.0-57.0 47.9 0.0359 55.1 9.4 


4.5 59 


2.0 


2.0 


1.5 


4.0 


12.0 


-2.9 


0.3 


-2.5 


0.1 


0.0 


-6.4 


57.0-80.0 68.3 0.0284 62.0 10.7 


5.3 57 


2.0 


2.0 


1.5 


4.0 


12.0 


-2.9 


0.3 


2.5 


0.1 


0.0 


-6.4 



Table 4: Elastic and proton-dissociative photoproduction cross sections a ((Ww) derived from 
the high- and low-energy data sets as a function of the photon proton centre-of-mass energy W 1P 
for the processes ep — \ e J/ip Y, where Y denotes either a proton p or a proton-dissociative sys- 
tem of mass m p < M Y < 10 GeV. These cross sections are obtained after the combination of 
the cross sections from the /i + fi~ and e + e~ decay channels and for the phase space as defined in 



table Hi (W?p) indicates the bin centres [55| and $^ is the transverse polarised photon flux per 
bin. At ot and A comb denote the total and the combined statistical and channel- specific individ- 
ual uncertainties, as obtained from the data combination, respectively. The global correlation 
coefficients p GC are also shown. The full covariance matrix can be found in [51]. The remain- 
ing columns list the bin-to-bin correlated systematic uncertainties corresponding to a + la shift 
due to the correlated tracking uncertainty (5™ corr , the correlated triggering uncertainty (55f ,corr , 



sys 
r2S 



the uncertainty from ip(2S) contributions S 2 y S , the integrated luminosities of the high- and low- 
energy data sets 5£**, Sfy , the tagging uncertainties in the LAr 5^ rl °, the plug <& UG and the 



FTS <5f y T s s , due to the modelling of the MC <^ s CModel , the Q 2 dependance <5g, and the cut on the 
empty calorimeter S^S. 
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Table 5: Elastic and proton-dissociative photoproduction cross sections derived from the high- 
energy data sets as a function of the squared four-momentum transfer at the proton vertex t, for 
the processes ep — > eJ/ip Y, where Y denotes either a proton p or a proton-dissociative system 
of mass m p < M Y < 10 GeV. These cross sections are obtained after the combination of the 
cross sections from the /i + /U~ and e + e~ decay channels and for the phase space as defined in 
table fl] (|t| bc ) indicates the bin centres p5J. The transverse polarised photon flux $^ for the 
given phase space range is 0.0953. See caption of table[4|for more dietails. 
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Table 6: Elastic and proton-dissociative photoproduction cross sections of the low-energy data 
sets as a function of the squared four-momentum transfer at the proton vertex t, for the processes 
ep — > J/ip Y, where Y denotes either a proton p or a proton-dissociative system of mass My > 
m p . These cross sections are obtained after the combination of the cross sections from the fi + fi~ 
and e + e~ decay channels and for the phase space as defined in tableHl (|t| bc ) indicates the bin 



centres [55 1. The transverse polarised photon flux $i, for the given phase space range is 0.1108. 



See caption of table pw for more dietails. 
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Figure 2: Di-lepton invariant mass distributions for the high- and low-energy data sets in the 
J/ip — > n + n~ decay channel, figures a) and b), respectively, and for the J/ip — > e + e~ decay 
channel, figures c) and d), respectively. For the muon sample the fits used to reconstruct the 
number of J/ip mesons are shown as well. For the electron sample the simulation of the QED 
background ep — > eX e + e~ is given by the shaded region and the J/ip signal and sideband 
normalisation regions are indicated. 



24 



a) 



b) 



o 

CD 



.9 0.8 



CD 

'en 

.CO 



c) 



— Elastic 

Proton dissociation 




40 50 60 70 80 90 100 110 
W m rec [GeV] 



o 

c 

CD 

"cd 

S 0.6 

h 0.4 

0.2 





1 



g 
'o 
2 0.8 

I 0.6 

.CO 



0.4 
0.2 





+ H1 data 
MC simulation 



40 50 60 70 80 90 100 110 



W 



YP, rec 



[GeV] 



— Elastic 

Proton dissociation 







d) 



g 
o 
2 0.8 

I 0.6 

.CO 



0.4 

0.2 





4 6 8 

-t rec [GeV 2 ] 



I H1 data 
MC simulation 



BV 







4 



6 8 

-t rec [GeV 2 ] 



Figure 3: Tagging efficiencies as functions of (a) W 1Jhrec and (b) — t rec as obtained from the 
simulations of elastic and proton-dissociative J/ip production. Tagging fractions as functions 
of (c) W 7P>r ec an d (d) —t rec , as obtained from the e + e~ data set in the invariant mass window 
m ee = 2.3 — 3.3 GeV. The data set contains elastic and proton-dissociative J/ip decays, as 
well as ep — > eX e + e~ events. It is compared to the simulation based on the event generators 
DIFFVM and GRAPE. The data (simulations) are shown by points (shaded bands). The vertical 
spread of the bands represents the uncertainty due to the tagging in the simulation. 
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Figure 4: Observed distributions as functions of W rPtrec and — t rec restricted in m u to the 
J/V signal region. The muon sample is shown in a) and b), the electron sample is shown 
in c) and d). The data, shown by the points, are compared to the simulation of elastic and 
proton-dissociative J/ip production. Also shown is the contribution from ip(2S) events and, for 
the electron sample only, the QED background. For the muon sample, background is subtracted 
from the data using a sideband method. 
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Figure 5: Combined elastic and proton-dissociative cross sections as a function of W 1P (circles) 
compared to the input data from J/tp — > e + e~ (triangles) and J/ip — > /i + /i~ (squares) of the 
HE and LE data sets. The error bars of the input data indicate the uncertainty composed of the 
statistical errors (inner error bars) and statistical errors combined with all individual systematic 
uncertainties (full error bars). The error bars of the combined data points reflect the uncertainty 
after the combination. The combined data points are drawn at their bin centres. The electron 
and muon data points are shifted in W 1P for better visibility. 
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Figure 6: Differential J/ip photoproduction cross sections da/dt as a function of the negative 
squared four-momentum transfer at the proton vertex, —t, as obtained in the high-energy data 
set for the (a) elastic regime and the (b) proton-dissociative regime and as obtained for the low- 
energy data set shown in (c) and (d). The error bars represent the total errors. Also shown by 
the curves is a simultaneous fit to this measurement and [|24J of the form da/dt = N e ie~ bel ^ for 
the elastic cross sections and da/dt = N pd (l + {b p d/n)\t\)~ n for the proton-dissociative cross 
sections. The fit uncertainty is represented by the spread of the shaded bands. 
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Figure 7: Proton-dissociative cross section as a function of —t (full circles) compared to a 
previous measurement at high |t| p4| (triangles) interpolated to match the W 1P , Q 2 and My 
ranges of the current measurement. The curve represents a simultaneous fit to both data sets, 
the spread of the shaded band its uncertainty. 
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Figure 8: J/ift photoproduction cross sections as a function of the photon proton centre-of-mass 
energy W JP for (a) the elastic and (b) the proton-dissociative regime. The data from the high- 
energy data set are shown by circles, the data from the low-energy data set as squares. The error 
bars represent the total errors. Shown by the curves is the simultaneous fit to the data from this 
measurement and [4|, see figure|9j The fit uncertainty is represented by the shaded bands. In (c) 
the ratio of the proton-dissociative to elastic J/ip photoproduction cross section is shown. The 
data are presented as full circles and the vertical bars indicate the total uncertainties, including 
normalisation uncertainties. The inner error bars represent the bin-to-bin uncorrelated errors, 
determined in an approximative procedure. The curve is the ratio of the fits shown in (a) and 
(b). The shaded band indicates the uncertainty on the ratio obtained from the fit uncertainties. 
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Figure 9: Elastic cross sections as a function of W 1P from this measurement compared to previ- 
ous measurements at HERA |4[[5j. The shaded band represents a fit to the present data and [J4j 
together with its uncertainties. 
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Figure 10: Compilation of elastic J/ip production cross section measurements including this 
measurement, previous HERA results [|4}[5j, results from fixed target experiments [26 27 1 and 



from LHCb [53 1. Also presented is the fit to the HI data only, indicated by the curve. The fit is 
extrapolated in W 1P from the range of the input data to higher values, as shown by the dashed 
curve. The shaded band indicates the uncertainty on the fit. 
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Figure 1 1 : Compilation of elastic J/ip production cross section measurements including this 



measurement, previous HERA results [|4}[5j, results from fixed target experiments [26 27 1 and 
from LHCb |j33). Also presented are QCD fits from [[14) to the previous HERA data J4|[2T|-|23 1 
to determine a gluon density at leading-order and next-to-leading order, indicated by the curves. 
The fits are extrapolated in W 1P from the range of the input data to higher values, as shown by 
the dashed curves. The shaded bands indicate the fit uncertainties. 
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